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I. INTRODUCTION 



Despite that no evidence of new physics signal has been observed at the colliders yet, 
the observations from cosmology reveal that more than 20% of the whole Universe is made 
up by the so-called dark matter (DM) [1]. But the Standard Model (SM) does not provide 
any candidate for the dark matter, so discovering dark matter would be an undoubtable 
evidence of new physics beyond the SM. Especially, a series of cosmic ray and gamma- 



ray observations from INTEGRA 
attracted extensively attention [7J, |8 
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27j, |28(, which indicate that: 



INTEGRAL detected a gamma-ray pick at 511 keV from the center of the galaxy, 
which can be explained with e + e~ annihilation there. 

ATIC observed-electron positron excess from 300 GeV to about 800GeV, while the 
observation from PAMELA shows that there is a positron excess at 10 — 100 GeV. 

Even though the bump-like structure observed by ATIC is not confirmed by Fermi 
LAT, an excess between 200GeV and ITeV still remains in the e + + e~ spectrum. 

DAMA/LIBRA and DAMA/Nal experiments reported positive results from the direct 
detection of dark matter 29]. 



Apparently, it is very exciting to interpret the anomalous signals as being induced by the 
dark matter. Toward this end, key questions that need to be answered are the following: 

• Why is there no hadronic anomaly detected in the cosmic ray experiments? 

• Where does the large boost factor (~ 10 2 — 10 4 ), which is necessary to interpret the 
positron anomaly as the dark matter signal, originate from? 

• What generates the pick near the electron mass pole observed by INTEGRAL? 

• Why other direct search experiments do not observe any signal of dark matter? 

All those mysterious questions might be related to different aspects of dark matter. Many 
dark matter models have been proposed to address on one or more questions listed above, 
but we expect there is a dark matter model which can explain all above experiments very 



Standard Model 
SU(2) L x U{1) Y 
W B $ 



FIG. 1: Sketch of our non-Abelian dark matter model 

well. In this work we propose a dark matter model in which the dark sector is gauged 
under a new SU (2) group. The dark sector consists of £77(2) dark gauge fields, two triplets 
dark Higgs fields, and two dark fermion doublets (dark matter candidates in this model). 
The SM sector interacts with the dark sector through kinetic and mass mixing operators. 
Using two component dark matters, the model could explain both PAMELA and Fermi LAT 
data very well, while satisfies constraints from both the DM relic density and SM precision 
observables. 

The paper is organized as follows. In Sec. II, we present the field content and the La- 
grangian of our model. In Sec. Ill, we discuss the astrophysics and cosmology observations. 
In Sec. IV, we show some constraints from the SM precision observables. A global fit to the 
experimental data is given in Sec. V. In Sec. VI, we analysis the signals of this model at 
the LHC. Discussions and conclusions are given in Sec. VII. The anomalous dimension and 
the evolution of the relevant operator, and the decay properties of some new particles are 
shown in Appendix A and B, respectively. 
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II. THE DARK MATTER MODEL 



In this section, we will propose a specific dark matter model. To answer the questions 
mentioned above, we adopt some ideas in references , |30|. Now we present a non- 



Abelian dark matter model (NADM) in details. As shown in Fig. [TJ our model consists of 
two sectors: one is the usual SM sector, the other is the dark sector (DS). The SM sector is 
gauged under the usual SU(3) C x SU{2) L x U(l)y symmetry while the dark sector under a 
new symmetry SU(2)d (named as dark symmetry throughout this work). Particles in the 
dark sector, including the dark matter candidates, neither carry any SM charge nor interact 
with SM fields directly. Similarly, all the SM particles are invariant under the dark gauge 
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group transformation. In other words, both sectors are invisible to each other, but they 
can communicate via (i) the mixing between the SM Higgs boson and dark scalar and (ii) 
high dimension operators induced by unknown ultraviolet (UV) complete theory at a much 
higher energy scale A. The Lagrangian of our models is 

£ = ^SM + ^DS + £mix, (1) 

where £sm (£ds) denotes the Lagrangian of the SM (dark) sector, and £ m j x represents 
the interaction between the SM and dark sectors. Here the dark matter candidates will 
annihilate into the dark gauge bosons and the dark scalar bosons, which eventually decay 
into the SM particles via mixing operators. 



A. Dark sector 



The dark matter candidates in our NADM are Dirac fermion doublets ipi and ip2- The 



most economic way to obtain anomaly free is to use the vector-like dark fermion. Two 
or more dark fermion doublets are needed in order to fulfill both XDM 7J] and iDM 30] 
scenarios, because both scenarios require different mass splits between the ground state and 
excited state of dark matter candidate. For example, the mass split is about 1 MeV in 
the XDM scenario but about 100 keV in the iDM scenario. In earlier works, some authors 
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321 ]. which does not need two 



suggest the dark matter fermion is dark gauge triplet 
generation DM fermions. 

It is worth emphasizing the possibility of multiple dark matter candidates. Many dark- 
matter candidates have been suggested in various models beyond the SM, but a nearly 
universal implicit assumption is that only one such candidate is needed and its properties 
are constrained accordingly. Of course, no fundamental principle requires there is only 
one dark matter candidate, and the possibility of multipartite dark matter should not be 



ignored [33, 
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381 ] . Moreover, as to be shown later, more than one generation 



dark matters can fit all experimental data much better in our NADM. 

The dark gauge symmetry could be broken in many ways, and one popular choice is with 
the help of scalar field as the Higgs boson doublet in the SM. However, as pointed out in 



Ref. 



391 ]. the iDM scenario cannot be realized due to a custodial symmetry when only one 



dark scalar doublet is added. Rather than adding two complex scalar fields, we use two real 
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scalar triplets (named as dark Higgs boson Hi) to break the custodial symmetry. 
The Lagrangian of the dark sector is 

3=1,2 

+tjj 1 (ip - mi) Vi + ^2 (ip - m 2 ) ^2 + E #7*^i%^i> (2) 

where the dark gauge field tensor J 7 ^, the dark Higgs field covariant derivative D^, and the 
fermion covariant derivative p are defined by: 

K> = d»Z« - d v Z; - g D e abc Z^Zt (3) 
D,H] = d,H] - g D e abc Z b d ^ (4) 
Pti = rfd^j + ign^Z^Y (5) 

Here, go is dark gauge coupling constant, Vh is the dark scalar potential, j3jk is the Yukawa 
coupling between dark Higgs boson and dark fermions, and rrii is the intrinsic DM fermion 
mass for ipi. 

We indicate that after symmetry breaking, both dark gauge bosons Z and dark fermions 
ipi obtain masses from the dark scalar vacuum condensation, and mz is proportional to 
gnVD while mj to PijVo- As to be discussed below, a large mass hierarchy between Z and 
ipi, say rriz <C rrii, is required to fit the PAMELA data. Hence, either a very light Z or 
a much heavy ipi is preferred. Unfortunately, the former requires a very small dark gauge 
coupling g£> which will over produce DM relic abundance, and the latter demands a huge 
Yukawa coupling which will destroy the vacuum stability. In order to avoid these problems, 
we explicitly keep the intrinsic fermion mass term m 8 in £ds- Such large masses could be 
generated by other exotic heavy fields which decouple at the scale much higher than the 
electroweak scale. The Yukawa interaction fijkipj'Hkipj would generates mass split between 
the two components inside one fermion doublet. 

Last, we comment on the dark scalar potential. When two and more scalar multiplets 
present, e.g. as in two Higgs doublet model, the scalar potential is rather complicated. 
Instead of discussing a general scalar potential, we will choose the scalar potential in Eq. [2] 
as a simple form, 

V H (Wi, H 2 ) = \ X (ll ( n i n i ~ ^/ 2 ) 2 " 2 W^) 2 j , (6) 
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FIG. 2: Representative Feynman diagrams which can induce the dimension-five operator: (a) heavy 
boson loop, (b) and (c) heavy fermion loop. 

which preserves a stable vacuum and breaks the custodial symmetry simultaneously. Here 
Vd denotes the vacuum expectation value (VEV) of dark Higgs bosons. 

B. Interaction between the SM and dark sectors 

We first consider that the dark sector interacts with the SM sector through a dimension-5 
operator 



where B^ u is the field strength tensor of the gauge boson associated with the SM U(l)y 
group and A is the cutoff scale of new physics (NP). Such an operator could be induced by 
integrating out new heavy particles in a renormalizable theory as shown in Fig. [21 which 
could generate the dimension-5 operator when new heavy particles decouple. The evolution 
behavior of this operator is determined by its anomalous dimension (see Appendix A). 

Another interaction between the SM and dark sectors is via the mixing between the dark 
scalar and the SM Higgs boson, such as 



where $ denotes the SM Higgs doublet and is the mixing parameter. Note that the scalar 
mixing is crucial to prevent overproducing DM relic abundance in our model. For example, 
without the scalar mixing, there exists a stable dark scalar after symmetry breaking, which 
could also be dark matter candidate as long as it is lighter than dark fermions. Even 



£ mix D -wtB^ (aiW? + a 2 H a 2 ) Z] 
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though the dark scalar pair could annihilate into the SM particles through the process 
TiTi. — > 77 —>■ SM particles, the annihilation cross section is too small to produce correct 
amount of relic abundance. In fact, it produces too much relic abundance which would 
overdose the Universe. Such a problem is solved by the scalar mixing with the SM Higgs 
boson, which enables the lightest dark scalar decaying into SM fermion. 
Next, we discuss some theoretical constraints on e^: 

• Vacuum stability : The existence of the lower bound of the scalar potential requires 

enVsMV d /fn H < M H , 

where t>sM and m# are the vacuum exception value and the mass of SM Higgs boson 
(if), respectively. 



Naturalness : As pointed out in Ref. 40J, the naturalness condition requires that 
e n v sM * s n °t muc h larger than v%, i.e. 



which leads to the following bounds: (1) en ~ 10 5 for Vd ~ 1 GeV and (2) ~ 10 7 
for v d ~ 0.1 GeV. 



C. Symmetry breaking and mass spectrum of dark particles 

In our model the dark SU(2)d gauge symmetry is broken spontaneously when the two 
dark scalar triplets develop non-zero vacuum expectation values (VEV). The dark scalar 
potential achieves its minimum value at 

f \ 





\ vd J 



and (H 2 ) = — 



Vd 

W 



(9) 



where, for simplicity, we choose same VEV (v^) for both scalar triplets. Therefore, the 
non-Abelian SU(2)d symmetry is totally broken by the non-zero Vd- 

After dark gauge symmetry breaking, all the three dark gauge bosons become massive: 



m Z2 = m Z3 = M z = g D v d , and m Zl = \/2M z . 



(10) 



The dark gauge boson mass will be modified slightly when they mix with SM gauge bosons 
through the operator as shown in Eq. [7J Without loosing generalization, we choose a\ — 
in Eq. [7J As a result, only Z 2 mixes with SM gauge boson. Note that a non-zero a\ merely 
changes the definition of mass eigenstates, and there is still only one dark gauge boson 
mixing with SM gauge boson directly. Then, the kinematic mixing between U(1)y gauge 
boson and Z 2 generated by this operator is 
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„B '^Z 2 , 
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where e = a 2 v d /y/2A, and = d^Z v 2 - d u Z%. It leads to non-diagonal elements in the 



kinetic energy term of Lagrangian in the basis v£ 



(Z^B^Wz,) 



K 
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where the gauge field with a caret symbol is understood as the current eigenstate. Further- 
more, the mass matrix takes the following form 
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We denote the mass eigenstates of the three gauge bosons by = (Z 2fl , A^, Z^). A simul- 
taneous diagonalization of both the kinetic energy term and the mass matrix gives (to order 

e) 



B» = A"c w - Z^sw + 



m z s w 
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where s w = sin9 w = g Y /Vg 2 + 9 Y , c w = cos9 w = # 2 /V#2 
corresponding mass eigenstates are given by (to order of e 2 ) 



g Y . And the masses of the 
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where rhz(=g2VsM /2cw) is the SM Z-boson mass. The gauge fields in the mass eigenstate 
basis are (to order of e) 



A» = c w B» + s w W£ - eZ%c w , 



Z» = c w W% - s w B^ 



_M\s w 
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m z ivi z 



(16) 



Two dark scalar triplets carry six degrees of freedom. While three of them are eaten by 
dark gauge boson after symmetry breaking, the rest three remain as three physical dark 
scalar fields. In the unitary gauge, the two scalar triplets can be written as 

/ n \ / „ \ 
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h 2 + v d /V2 
h 3 /V2 



(17) 



where hi, h 2 and h 3 are physical dark scalars with degenerate mass = Xdv^ = <\dM| 

In addition the mixing term in Eq. [8] can induce the following mass matrix in the basis (Hq, 

h 3 ) 



\ ( Ho h 
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Thus, the mass eigenstates and their masses are given by 
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l h ivi n m H 

The dark Higgs boson h 3 is no longer stable as it can decay to SM fermion pair through the 
above mixing. 

III. ASTROPHYSICS AND COSMOLOGY OBSERVATIONS 
A. Dark matter relic density 



M 2 H 



(19) 



In our model the dominant contribution to dark matter annihilation comes from the 
channels involving dark gauge boson and dark Higgs (also with the lighter generation dark 
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fermion for the heavier generation freezing out) in the final state, which is of order a 2 D . 
Those annihilation products decay into the SM particles eventually. Below we calculate the 



relic abundance in our model following Ref. |41|. Furthermore, coannihilation effects [42J are 
also considered in our calculation, as the masses of the two component fields of DM fermion 
doublet are nearly degenerated. 

As to be shown later, in order to fit both PAMELA and Fermi LAT data in our model, it is 
necessary to have a large mass splitting between ipi and ip2- Therefore, the DM annihilation 
involves the dark gauge interaction at different energy scales. Running effect of dark gauge 
coupling has to be considered in the numerical evaluation. The running behavior of dark 
gauge coupling is governed by renormalization group equation (RGE). The one-loop level 



beta function for a general non-Abelian gauge theory is 



43| 



^> = i&(-Y" + ^ + H' (20) 

In our model, N = 2, C2& = 2n b and C2/ = , where n b = 2 is the number of real Higgs 
triplet, n/ is the number of fermion doublet. For the energy scale /1 > 777.2, n f — 2, for 
mi < fi < 777.2, Uf — 1, and for fi < mi, rif = 0. Within this theoretical framework, the 
dark gauge interaction is asymptotic free, i.e. becoming stronger at low energy but weaker 
at high energy. In this work we choose the following parameters as benchmark: 

mi = 510 GeV , m 2 = 1300 GeV, 

M z = 0.1 GeV , M n = 0.02 GeV, (21) 

to explore the impact of RG running effect on the relic abundance. 

Figure [3]^a) displays the relic density versus dark gauge boson coupling at the scale 
\i = Mz- Note that, with running effects, the dark gauge coupling is evaluated at the scale 
fi = 2m i in the numerical calculation of DM annihilation cross sections. It is well known that 
the relic density of DM is approximately given by flh 2 ~ (0.1 pb)/ (ay) where (av) is the 
thermally averaged product of the DM annihilation cross section with its velocity. Five-year 
WMAP data indicates (av) ~ pb. Since the final state particles in the DM annihilation 
are almost massless, the annihilation cross section can be written as following, based on the 
dimension counting, 

<TiV oc (22) 
mj 
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When the DM candidate masses are given, there exists an unique dark gauge coupling to 
produce the correct amount of relic abundance. For example, without running effect, the 
correct relic density can be produced only in the vicinity of old ~ 0.1; see the black curve 
in Fig. [3](a). With running effects, the gauge coupling at the scale \x = 2m i is less than its 
value at the scale /i = M z , therefore, large coupling strength is needed to induce efficient 
DM annihilation; see the red curve in Fig. [3]^a). Owing to the RG running effect, the favored 
coupling region, which is consistent with the WMAP data at la level (green band), become 
much broader (e.g. 0.45 < ap < 0.67) for the given benchmark parameters. 

In Fig. E](b), we display the comparison of two dark matters (red) and one dark matter 
(black). In the case of one generation dark matter, we fix its mass be 1300 GeV because 
the Fermi LAT data prefers a very heavy DM candidate. It is clear that the model cannot 
produce a correct DM relic density when only one heavy DM candidate presents, even for 
an unreasonably large coupling such as ap ~ 5. It is due to the heavy mass suppression in 
Eq. [22] which yields a small DM annihilation rate, cf. the black curve. Introducing another 
generation dark matter efficiently reduces relic density to satisfy the WMAP data, cf. the 
red curve. 



B. INTEGRAL and XDM 



As pointed out in Ref. [7j, if dark matter has an excited state and the mass split is 
about 2m e , decay of the excited state of dark matter may produce much e + e~ as a source 
of the gamma-ray line in the galactic center at 510.954 ± 0.075 keV confirmed by INTE- 
GRAL/SPI [44] • Following the proposal in Ref. [7|, we find that in our model the dark gauge 
boson mass Mz should be less than 0.2 GeV in order to give a significant contribution to 
the 511 keV signal. As a result, such a light dark gauge boson cannot decay into quarks, 
muon and tau due to the kinematics. 



C. PAMELA and Fermi LAT 

To fit the data of PAMELA and Fermi LAT, we calculate the cosmic ray e + e~ flux, using 



the background flux formula 45| and the formulas of the dark-matter- induced flux [46!]. In 



Appendix B we summarize the decay properties of dark gauge boson and dark scalars used 
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FIG. 3: (a) Relic abundance versus dark gauge coupling evaluated at the scale Mz, with and 
without running effect of dark gauge boson. The annihilation cross sections are calculated with 
the dark gauge coupling at the scale /j = 2m«. The parameters are chosen to be m\ = 510 GeV, 
m 2 = 1300 GeV, M z = 0.1 GeV and M H = 0.02 GeV. The dark green band display the cold dark 
matter density from Five- Year WMAP data, £lh 2 = 0.1099 ± 0.0062 [1]. (b) Comparison of two 
generation DMs and one generation DM with the choice of above parameters, except that we fix 
the DM mass be 1300 GeV in the one dark matter case. 
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in our numerical calculation. The large Sommerfeld enhancement factor, which is give n by 
the non-Abelian dark gauge interaction, can be calculated numerically following Ref. 47]. 
The results of fitting PAMELA and Fermi LAT is shown below. 



D. Direct detection of dark matter 



Current direct detection experiments of dark matter 



like 



23, 



48 



49 



50] give negative result 



except DAMA 29J]. In our model, the interaction between dark matters and nucleon is iDM 



30], so the direct detection can not give a strong limit to e. For the generation whose 



mass split is greater than 2m e required by XDM, the corresponding critical velocity [3( 
will be larger than the escape velocity of dark matter in the galaxy, which means that it 
can not be detected in the direct detection experiments. For another generation, the mass 
split is still a free parameter which can be use to explain DAMA and other direct detection 
experiments data. 



IV. CONSTRAINTS FROM THE SM PRECISION OBSERVABLES 

Due to its mixing with the SM gauge bosons, the dark gauge boson inevitably modifies 
various low energy precision observables, such as the oblique parameters , muon anoma- 
lous magnetic momentum, and Z° boson decay, etc. In this section we explore constraints 
from those precision measurements. 

a. The Electroweak Oblique Parameters The oblique parameters S, T and U can be 



calculated with formulas in Ref. 



521 ] . In our model, to the order of e 



S 



G F (ml-Mp 1 



T 



Zi 
2 



V2c w G F (m 2 z -M 2 z Y 
U = 0, (23) 

where Gp is the Fermi constant. For Mz ~ 0.1 GeV, from above parameters we can get 
e < 0.023 at lex level. 
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b. Muon Anomalous Magnetic Moment The additional contribution to the muon 
anomalous magnetic moment (g — 2) is mainly from Z2, which can be calculated by 53] 



F 2 (0) = MmJ fdzil-z) I 
m Zo Jo J 



with 



(24) 



h x {z) = -4z(z + 3), 

fc 2 (*) =4z(l-z), 

2 

6 2 = z + ^(l-z) 



Z2 



and 



9l 



e 2 e 2 
\Qc 2 w 
e 2 e 2 M| 



4:C 2 W + 



Ml 
M 2 



(25) 
(26) 
(27) 



(28) 



where m M is the muon mass. For e < 7 x 10 4 , the contribution from Z<± to 



— 2)„ is less 

than 1(T 10 . Here we choose M z ~ O.lGeV and M n ~ 0.02GeV. 

The contribution from /i 3 is much less than iJ 2 . Even for ~ 400 the contribution from 
h 3 to (g 

show that the constraint from 

c. Decay width of the SM Z° Boson In our model there are some new decay channels 
of the SM Z°: Z° -> Z ± Z 3 , Z° -> Z 2 /i 1; Z° -> Z 2 /i 2 , Z° -> Z 2 /i3, Z° -> 2^3, Z° -> 7 /i 2 . 



2) is less than 10~ 10 when M z ~ O.lGeV and M H ~ 0.02GeV. These results 

— 2) is not strong. 



Ignoring the higher order terms of M^/m 2 ^ and M\jm 2 Zl to the order of e 2 , the non-zero 



decay widths are 



1 z r z 3 



1 Z 2 h 2 



2rf r 

Z3I13 



Qpe 2 s 2 y m| 
12M| ' 



(29) 



1 Zlhx 



a D e 2 s 2 v m z 
24 ' 



(30) 



From above results and the experimental data of the Z° decay width 54], at 1 a level, e 
should be less than 4.5 x 10~ 5 for M z ~ O.lGeV, a D (M z ) = 0.57 and M H ~ 0.02GeV. 
Our calculation shows that the evolution of the operator in Eq. [7] is negligible small (see 
Appendix A). 
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V. EXPLANATION OF THE PAMELA AND FERMI LAT OBSERVATIONS 



In this section, we present a global fit of the PAMELA, Fermi LAT, DM relic density, the 



S and T parameters, and Z° boson decay width, using least y 2 analysis 55|, |56j . As widely 
used in the literatures, a parameter 6 fc is introduced to rescale the cosmic ray background 
here. In general, the total distribution of dark matter in the galaxy can be determined 
by N-body simulations, but for a model of two generation dark matters, the fraction of 
each generation in the galaxy is unknown. In our model, since the masses of two dark 
matters are different, the ratio r ga i axy of their total masses in the galaxy will be different 
from rf r ee ze = ^heavy /^Hght > where r^heavy (^Hght) is the relic density of the heavier (lighter) 
dark matter. So we take r galaxy as a free parameter. 

In our fitting we require the theoretical values of both the SM relevant precision observ- 



ables and the relic density fall in la region o: 



and Fermi LAT data using NFW profile [57] and "MED" propagation scheme 



the experimental data. We fit the 



'AMELA 



46|. The 



best-fit solution is a D (M z ) = 0.57, M z = 0.1 GeV, M H = 0.02 GeV, m x = 510 GeV, 
m 2 = 1300 GeV, e = 10~ 5 , bk = 0.66 and r ga i axy = 6.1, resulting in \ 2 — 11 P er degree of 
freedom. 

We illustrate, in Fig. BJa), the positron fraction as seen at Earth after propagating effects 
are included in the "MED" propagation scheme 4(|. We note that our best fit is quite 
good, with the Sommerfeld enhancement boost factor of ~ 5 x 10 3 in our case. The kink 
at the high energy region is due to the overlap of cosmic positrons from two component 
DM annihilations. Figure IU(b) shows that our model prediction also fits the Fermi LAT 
very well. With enough precision, the kink feature of two component dark matters may be 
explored in future experiments. 



VI. LHC PHENOMENOLOGY 



The typical signals of our model at the LHC are the lepton jets, which are boosted groups 
of n > 2 leptons with small angle separations and GeV scale invariant masses as shown in 
Ref. [66J. Especially, for the process p + p — > Z° — > Z 2 + h 2 — > e + e~ + $t, the signal is just 
one lepton jet originated from the Z 2 decay plus $t from the long lived h 2 . See Appendix IE1 
for detailed calculation of the dark gauge boson and dark scalar decays. The intrinsic SM 
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backgrounds originate from the W + W~ and ZZ/Z^f* pair production. Only the W — > ez/ e 
decay mode is considered in the WW background in order to mimic the signal signature. 
Similarly one of the Z boson has to decay into neutrino pair and the other Z boson or off- 



671 ]. we examine 



shell photon decays into electron-positron pair. With the help of CalcHEP 
the kinematics distributions of both signal and background. 

Kinematics of the signal process is very distinctive from those of background processes: 

• Signal: The two charged leptons in the signal are forced to move collaterally due to 
the large boost received from Z 2 . Even though it is challenging to measure the mo- 
mentum of each charged lepton, one can measure the sum of the energy and transverse 
momentum (pt) of the two charged lepton system. The distribution of p^ e as well as 
$t peaks around half of the Z° boson mass; see the red-dotted curve in Fig. [5j Similar 
to the W — ► Iv decay, the transverse mass of the lepton-jet and missing transverse 
momentum J£t also exhibit a sharp Jacobin peak around m z o. On the other hand, 
the invariant mass of the leptons is around GeV. Last, the Z 2 boson is long lived with 
e = 10~ 5 . It, when produced, will propagate about 30 centimeter inside the detector, 
which provides an unique collider signature. 

• WW background: The two charged leptons in the WW background neither move 
parallel to each other nor exhibit GeV invariant mass peak. 

• ZZ background: Since the dominant contribution is from the on-shell Z boson pro- 
duction, the invariant mass of two charged leptons is close to m^o. The distribution 
of p^ e , i.e. peaks around 30 GeV. 

• Z--/* background: It provides very similar collider signature as the signal. Originated 
from a off-shell photon decay, two charged leptons move in parallel and have a small 
invariant mass. The transverse mass of the lepton system and I£t does not peak 
around m z o. Unlike the signal, the off-shell photon decays promptly and do not travel 
a long distance inside the detector. 

The SM backgrounds, mainly from WW production, overwhelm the signal. Making use of 
the kinematics difference mentioned above, one can impose optimal cuts to suppress the SM 
background. For illustration, we impose the following two simple cuts: 

m e+e - < 10 GeV and cos# e+e - > 0.95, (31) 
17 
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FIG. 5: Differential distribution of transverse momentum of the lepton jet system p 1 ^ 6 with the 



cuts: m P +„- < lOGeV and cos i 



> 0.95. 



where m e + e - is the invariant mass of lepton-jets and 9 e + e - is the open polar angle between 
the electron and positron. Figure [5] shows the p^ e distributions of signal and background 
with the cuts in Eq. [3H which clearly shows that it is very promising to observe the dark 
gauge boson signal at the LHC However, one should bear in mind that the above two cuts 
serve for the purpose of suppressing the SM background processes. It is crucial to understand 
how well one can measure the energy and momentum of the lepton jets, which is beyond 
the scope of this paper and will be presented in the future work. 



VII. CONCLUSIONS 



Discovering dark matter would be an undoubted evidence of new physics beyond the 
Standard Model. Recently, many excesses were observed in the direct and indirect search of 



dark matter 



29j . In this work we proposed a non-Abelian dark matter model to 



explain all those observed excesses. The model consists of both SM sector and dark sector. 
The latter, gauged under a new SU(2) dark gauge symmetry, contains SU{2) dark gauge 
boson fields, two triplets dark Higgs fields and two generation dark fermion doublets (dark 
matter candidates). Rather than considering only one dark matter candidate, we introduced 
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two coexisting stable dark matters: one is around 500 GeV and the other is around 1300 GeV. 
Our study shows that two dark matter scenario naturally fits both PAMELA and Fermi LAT 
data simultaneously. Dark matter annihilation occurs in the dark sector and the remanent 
dark particles (dark gauge boson and dark scalars) decay into the SM particles eventually 
through the kinetic and mass mixing operators between the SM and dark sectors. The triplet 
dark scalars are introduced to break the dark gauge symmetry and also generate small mass 
splits between the two component fields of dark fermion doublets. Such small mass splits is 
the key to realize the iDM and XDM scenarios, which explain the INTEGRAL and DAMA 
results. Finally, we explore the interesting collider signature of the dark gauge boson and 
dark scalar production at the LHC. Our simulation analysis indicates that the signals of this 
model may be detectable at the LHC. 

In our model, since the Lagrangian of the dark sector has scalar mass terms, there is a 
dark sector fine tuning problem. In general, such a fine tuning problem can be solved by 
introducing supersymmetry in the dark sector. The supersymmetry extension of this model 
will be especially interesting and will be presented elsewhere. 
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APPENDIX A: ANOMALOUS DIMENSION OF THE OPERATOR 

In this appendix, we show the evolution of the dimension-five operator in Eq. [7J since 
the operator evolution can change the theoretical result significantly [68j. The evolution of 
the Wilson coefficient C{p) of the effective Hamiltonian ±C(ji)0(p) (O = B^Z^Uf, is 
determined by the evolution equation 

Hj-C = ^oC, (Al) 
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FIG. 6: The feynman diagrams which contribute the operator anomalous dimension. 

where 70 is the anomalous dimension. Calculating the Feynman diagrams in Fig. 6, we can 
obtain 



9l 



24tt 2 



(1-n/) 



(A2) 



where rif is the number of fermion doublet as defined in Sec. III. A. The solution of the 
evolution equation (Al) is 



C (ji) = exp 



9»Q>) lo (g 



D, 



dg 



D 



(7(A). 



lg D (A) P(9d) 

Thus, the Wilson coefficient at Mz scale becomes 

c r Mz \ 1 a D{M z )\ ( a D {m 2 )\ 



(A3) 



\ a D (mx) J \ a D (A) 



C(A). 



(A4) 



APPENDIX B: DECAY OF DARK GAUGE BOSONS AND DARK SCALARS 



In this appendix, we discuss the decay properties of dark gauge bosons and dark Higgs 
bosons. 
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1. decay of Z\ 



The main decay channels of Z\ are Z\ — > h\hi and Z\ — > /11/13. The total decay width is 
given by 

For M 2 = O.lGeV, M w = 0.02GeV and « D (M 2 ) = 0.57, we have T 2l ~ 0.006GeV. 



2. decay of Z 2 

The main decay channels of Z 2 are Z 2 — > 7/12 and Z 2 — > // (f — e~, v e , u T ). For 
7/2,2 channel, the decay width is 

r^ 2 - — T2 — [i-jq) ■ ( B2 ) 

The decay width of channel Z 2 — * e + e~ is 



_ ae 2 c 2 w M z ( 5m 2 e m\ \ Am 2 e a> . r] 



where a is the fine structure constant. The decay widths T Z2 _ Uii7i , (i = e, //, r) are negligible 
small. So the total decay width of Z 2 can be written as 



ae 2 c 2 w M z ( 5m 2 e m\ \ 4m 2 e a D e 2 c 2 w M z ( M 2 n V 
Tz > = 3 \ l + 2Ml~2Ml) V Ml + 12 V M| / ' 

For M 2 = O.lGeV, M H = 0.02GeV, e = 10" 5 and a D (M z ) = 0.57, we have Y z , 2 ~ 3 x 
10- 13 GeV. 



3. decay of Z 3 



The main decay channels of Z^ is Z 3 — > 7/13. The total decay width is 



_ ap ^M z / M 2 \ 3 
r * _ 24 V M|J ' (B5) 



For M z = O.lGeV, M H = 0.02GeV, and a D (M z ) = 0.57, we have T Z . A ~ 1.6 x lQ- 13 GeV. 
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4. decay of hi 



The main decay channels of hi are hi — > fff'f where /, /' = e~, v e , v T . For 
M z = O.lGeV, M H = 0.02GeV and a D (M z ) = 0.57, the numerical result of the decay 
width shows that it is of order O(10~ 35 GeV). 

5. decay of h 2 

The dominant decay channels of hi are h 2 — > 7// (f = e~ , v e , u T ). The numerical 
calculation shows that Fh 2 ^ 1Uii7i , (i = e, r) are negligible small. Neglecting the electron 
mass and the contribution from diagrams with the Z° boson in the internal line, the total 
decay width can be written as 



h2 ~ 72n 



Ml M% ( 4Mf\ / M\\ 2 ( Ml 

-17— g + 42 - 24—| + 61- — # 1 - — # In 1 n 



M 2 Z M^ \ Ml J V Ml J V M z 



(B6) 

For M z = O.lGeV, M z = 0.02GeV, e = 1(T 5 and a D (M z ) = 0.57, r (l2 ~4x 10~ 32 GeV. 
6. decay of J13 

h 3 decays into e + e~ . The decay width is 

e 2 n M n m\Ml f _ 4r^ 
32vr 2 a D \ M w 




r h 3 ^^2— li Z ( l - TJt) Xl 1 - T7T- ( BT ) 



For m H = 115GeV, M z = O.lGeV, M H = 0.02GeV, e H = 10~ 7 and a D (M z ) = 0.57, the 
decay width of h 3 is 1.7 x 10~ 35 GeV. 



[1] J. Dunkley et al. (WMAP), Astrophys. J. Suppl. 180, 306 (2009), 0803.0586. 

[2] A. W. Strong et al., Astron. Astrophys. 444, 495 (2005), astro-ph/0509290. 

[3] J. Chang et al., Nature 456, 362 (2008). 

[4] O. Adriani et al. (PAMELA), Nature 458, 607 (2009), 0810.4995. 

[5] O. Adriani et al., Phys. Rev. Lett. 102, 051101 (2009), 0810.4994. 

[6] A. A. Abdo et al. (The Fermi LAT), Phys. Rev. Lett. 102, 181101 (2009), 0905.0025. 



22 



[7] D. P. Finkbeiner and N. Weiner, Phys. Rev. D76, 083519 (2007), astro-ph/0702587. 
[8] L. Bergstrom, T. Bringmann, and J. Edsjo, Phys. Rev. D78, 103520 (2008), 0808.3725. 
[9] M. Cirelli and A. Strumia (2008), 0808.3867. 
[10] V. Barger, W. Y. Keung, D. Marfatia, and G. Shaughnessy, Phys. Lett. B672, 141 (2009), 
0809.0162. 

[11] J.-H. Huh, J. E. Kim, and B. Kyae, Phys. Rev. D79, 063529 (2009), 0809.2601. 

[12] N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer, and N. Weiner, Phys. Rev. D79, 015014 

(2009), 0810.0713. 
[13] P.-f. Yin et al., Phys. Rev. D79, 023512 (2009), 0811.0176. 

[14] K. Ishiwata, S. Matsumoto, and T. Moroi, Phys. Lett. B675, 446 (2009), 0811.0250. 
[15] A. Ibarra and D. Tran, JCAP 0902, 021 (2009), 0811.1555. 

[16] D. Hooper, A. Stebbins, and K. M. Zurek, Phys. Rev. D79, 103513 (2009), 0812.3202. 
[17] Q.-H. Cao, E. Ma, and G. Shaughnessy, Phys. Lett. B673, 152 (2009), 0901.1334. 
[18] P. Meade, M. Papucci, and T. Volansky (2009), 0901.2925. 
[19] D. Hooper and K. M. Zurek, Phys. Rev. D79, 103529 (2009), 0902.0593. 
[20] C. Cheung, J. T. Ruderman, L.-T. Wang, and I. Yavin, Phys. Rev. D80, 035008 (2009), 
0902.3246. 

[21] K. Cheung, P.-Y. Tseng, and T.-C. Yuan, Phys. Lett. B678, 293 (2009), 0902.4035. 
[22] X.-J. Bi, X.-G. He, and Q. Yuan, Phys. Lett. B678, 168 (2009), 0903.0122. 
[23] S.-L. Chen, R. N. Mohapatra, S. Nussinov, and Y. Zhang, Phys. Lett. B677, 311 (2009), 
0903.2562. 

[24] K. Ishiwata, S. Matsumoto, and T. Moroi (2009), 0903.3125. 
[25] A. A. El-Zant, S. Khalil, and H. Okada (2009), 0903.5083. 

[26] V. Barger, Y. Gao, W. Y. Keung, D. Marfatia, and G. Shaughnessy, Phys. Lett. B678, 283 
(2009), 0904.2001. 

[27] P. Meade, M. Papucci, A. Strumia, and T. Volansky (2009), 0905.0480. 

[28] X.-G. He (2009), 0908.2908, and references therein. 

[29] R. Bernabei et al. (DAMA), Eur. Phys. J. C56, 333 (2008), 0804.2741. 

[30] D. Tucker-Smith and N. Weiner, Phys. Rev. D64, 043502 (2001), hep-ph/0101138. 

[31] F. Chen, J. M. Cline, and A. R. Frey, Phys. Rev. D79, 063530 (2009), 0901.4327. 

[32] F. Chen, J. M. Cline, and A. R. Frey (2009), 0907.4746. 

23 



[33] E. Ma, Annales Fond. Broglie 31, 285 (2006), hep-ph/0607142. 

[34] T. Hur, H.-S. Lee, and S. Nasri, Phys. Rev. D77, 015008 (2008), 0710.2653. 

[35] Q.-H. Cao, E. Ma, J. Wudka, and C. P. Yuan (2007), 0711.3881. 

[36] H. Sung Cheon, S. K. Kang, and C. S. Kim (2008), 0807.0981. 

[37] M. Fairbairn and J. Zupan, JCAP 0907, 001 (2009), 0810.4147. 

[38] S. Profumo, K. Sigurdson, and L. Ubaldi (2009), 0907.4374. 

[39] M. Baumgart, C. Cheung, J. T. Ruderman, L.-T. Wang, and I. Yavin, JHEP 04, 014 (2009), 
0901.0283. 

[40] D. P. Finkbeiner, N. Padmanabhan, and N. Weiner, Phys. Rev. D78, 063530 (2008), 
0805.3531. 

[41] E. Kolb and M. Turner, Addison- Wesley P.C. (1990). 
[42] K. Griest and D. Seckel, Phys. Rev. D43, 3191 (1991). 
[43] D. J. Gross and F. Wilczek, Phys. Rev. D8, 3633 (1973). 

[44] E. Churazov, R. Sunyaev, S. Sazonov, M. Revnivtsev, and D. Varshalovich, Mon. Not. Roy. 

Astron. Soc. 357, 1377 (2005), astro-ph/0411351. 
[45] E. A. Baltz and J. Edsjo, Phys. Rev. D59, 023511 (1998), astro-ph/9808243. 
[46] M. Cirelli, R. Franceschini, and A. Strumia, Nucl. Phys. B800, 204 (2008), 0802.3378, and 

references therein. 
[47] R. Iengo (2009), 0903.0317. 

[48] G. J. Alner et al., Astropart. Phys. 28, 287 (2007), astro-ph/0701858. 
[49] J. Angle et al. (XENON), Phys. Rev. Lett. 100, 021303 (2008), 0706.0039. 
[50] Z. Ahmed et al. (CDMS), Phys. Rev. Lett. 102, 011301 (2009), 0802.3530. 
[51] M. E. Peskin and T. Takeuchi, Phys. Rev. D46, 381 (1992). 

[52] K. S. Babu, C. F. Kolda, and J. March-Russell, Phys. Rev. D57, 6788 (1998), hep-ph/9710441. 

[53] I. Bars and M. Yoshimura, Phys. Rev. D6, 374 (1972). 

[54] C. Amsler et al. (Particle Data Group), Phys. Lett. B667, 1 (2008). 

[55] D. Stump et al., Phys. Rev. D65, 014012 (2001), hep-ph/0101051. 

[56] C. S. Li, Z. Li, and C. P. Yuan, JHEP 06, 033 (2009), 0903.1798. 

[57] J. F. Navarro, C. S. Frenk, and S. D. M. White, Astrophys. J. 462, 563 (1996), astro- 
ph/9508025. 

[58] M. Aguilar et al. (AMS-01), Phys. Lett. B646, 145 (2007), astro-ph/0703154. 

24 



[59] R. L. Golden et al. (WIZARDS), Astrophys. J. 457, L103 (1996). 

[60] M. Boezio et al, Astrophys. J. 532, 653 (2000). 

[61] S. W. Barwick et al. (HEAT), Astrophys. J. 482, L191 (1997), astro-ph/9703192. 

[62] K. K. Tang, Astrophys. J. 278, 881 (1984). 

[63] M. A. DuVernois et al., Astrophys. J. 559, 296 (2001). 

[64] S. Torii et al., Astrophys. J. 559, 973 (2001). 

[65] F. Aharonian et al. (H.E.S.S.), Phys. Rev. Lett. 101, 261104 (2008), 0811.3894. 

[66] N. Arkani-Hamed and N. Weiner, JHEP 12, 104 (2008), 0810.0714. 

[67] A. Pukhov (2004), hep-ph/0412191. 

[68] E. Braaten, C.-S. Li, and T.-C. Yuan, Phys. Rev. Lett. 64, 1709 (1990). 



25 



